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We study the full set of planar Green's functions for a two-matrix model using the language 
of functions of non-commuting variables. Both the standard Schwinger-Dyson equations 
and equations determining connected Green's functions can be efficiently discussed and 
solved. This solution determines the master field for the model in the 'C-representation.' 
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A recent convergence of work by mathematicians [|T],0 and physicists has 
produced precise constructions of the 'master field' for general large N gauge and matrix 
field theories. The essential difficulty of the large number of degrees of freedom in higher 
dimensional large theories is dealt with by finding master fields which live in 'large' 
operator algebras such as the type IIi factor associated with a free group. 

When we say that the master field was 'constructed' in a higher dimensional theory, so 
far what we mean is that if we know all correlation functions for the theory, we can produce 
an operator representation of the master field. Thus to begin we should study models we 
can already solve, and see if any simplifications appear in the new language. The large 
number of degrees of freedom of higher dimensions (and the need to use non-trivial free 
algebras) appears in the large N limit of any integral over more than one matrix, and the 
simplest example is the two-matrix model defined by the large N limit of the integral 

VMi VM2 e~^^ [V{M^)+V{M2)-cM^M2] 

with Mi hermitian. We will take V to be the cubic polynomial V{M) = \M'^ — |M'^ 
below, but our considerations generalize easily to higher order polynomials. This model 
was first solved in using orthogonal polynomials. This technique uses in an essential 
way the special form of the action in and the ability to express the solution in terms 
of the invariant observables constructed from a single matrix. This generalizes to higher 
dimensions only in very special cases such as the Kazakov-Migdal model [jl2| and we would 
like to explore techniques which do not depend so fundamentally on a particular form for 
the action. 

One such technique is to solve the factor ized Schwinger-Dyson equations. At present, 
one can do this explicitly only when one can find a subset of these equations which closes on 
a few invariants, and such a reduction for ([^) was found in [lTO| , pT[] . We begin by discussing 



this and the recursive definition of higher correlation functions in the language of . The 
generating function of planar Green's functions is a function of non-commuting variables 

^{ui,U2)= ^ w{ui,U2) {^Trw{Mi,M2)). (2) 

words w 

The Green's functions and thus (j) have cyclic symmetry. 
The Schwinger-Dyson equations for then read 

did) — cd2(p — =4>ui4> 

(3) 

d2(j) — cdi4> — gd2(j) =(j)U2(p- 
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Derivatives acting on functions of free variables satisfy 

d 



dm 



Ujf{u) = diUjf{u) = 5ijf{u) (4) 



and this is just an alternate notation di = ai, Ui = a* for the free algebra of [|1|. We also 
define right derivatives 

/(w)wj^ = f{u)ujdi = 5ijf{u). (5) 
Acting on a function with cyclic symmetry (such as (f)) these satisfy 

d...f{u)% = d...d,f{u). (6) 

To get a truncated system of equations we work with the first terms in the series 
expansion 

^2<^(«l5'f^2)|„2^o = (/>o(wi) + U2(t)l{ui) + ul(t)2{ui) + . .., (7) 

consider the equations (^ and their derivatives at U2 = 0, and look for a minimal set of 
closed equations. Identifying terms of the form (|^) in (^ we get 

di(j)k - c(j)k+i - gdl4>k =(j)kUi4'o (8) 
- cdi(t)o - g(p2 =0. (9) 

From (P) a minimal closed set will be {^o, 4>ii4'2} and including (H) with k = and k = 1 
gives three equations. The operator di is then rewritten 
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n-l 



9ifM = [-^fiui)h = - I fM - ^ u\ dlfiO) I (10) 



i=0 



producing algebraic equations with boundary conditions 9{0(O) for i < 3. These can be 
combined into a cubic equation for (po{ui), to be made more explicit below. In general, we 
get an equation of degree degV. 

This is just a start on the problem of getting all correlation functions and the simplest 
description of these would be a single equation for (p which 'closes.' Not having very 
sophisticated techniques for working with functions of non-commuting variables, what we 
will mean by this is that we can directly truncate this single equation to an equation for 
^o(^^i) and then solve iteratively for cp in an expansion with terms containing U2^s in all 
possible ways. 
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We can find such an equation by repeating the strategy of taking derivatives of (||) 
and eliminating common factors in the result, but this time not setting U2 = 0. A useful 
set of equations derived from is 

d2di(j) - gd2df(j) - - {d24> - cdi(p - (j)U2(j)) = d2(l)Ui(j), (11) 

and equations derived by acting with the right derivative ^ 

d2(f) =- (di(p- gdf(p - (f)Ui(p) 
c 

d2di(P =^ {df(P - gdf(P - (Puid^P - (P) (12) 
c 



The final equation is obtained by substituting eqs.(^2l) into (|TT| 
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C , , C 



D(j) = (f)U24> 4'Ui(f) + di(f)Ui(f) + (j)Uidi(f) — gdi4'Ui(f) — g(f)Uidi(f) ~ 4'Ui(f)Ui(f). (13) 

where the operator D is given by 
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D^-l + gd,m + (9-- + -)di + (1 + c)df - 2gd! + g^df. 
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First let W2 = in (p!3|), to obtain a cubic equation for 0o(wi). Let R{z) — l/z(j)o{l/z), 
this algebraic equation is re-written as 

F(z, R) = R^ - f{z)R^ + g{z)R - h(z) = 0, 

f{z) =-- + 2z- 2gz\ 
9 

g{z) ^g'z' - 2gz'' + (1 + c)z' - {g + c/g - d'/g^)z + 1 - gd^m. ^^^^ 
h{z) =g^z^ + {g^d^m - 2g)z^ + {g^dlm - 2^^10(0) -l-c)z 
+ g'^dlm - 2gdlm + (1 + c)d^m +g-c/g + / g\ 

The 'initial' data 9"^(0) are not totally independent, as can be seen from eq.(|D. They 
are determined by requiring the appropriate analytic behavior of R{z), which for small g 
must have a single cut on the real axis. This equation was derived in (see also [1T0[]). 
In principle one can solve (|13D iteratively for (p in an expansion 

</>=5^<^(")(z.i,^.2), (15) 

n=0 
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where (j)^'^^ contains all terms with n instances of U2- Since there is only one term involving 
U2 explicitly, all terms containing cp^"^^ in the n-th recursion relation are linear in (p^'^^ and 
depend only on ^o- One can show that this differential operator, linear in (f)^'^\ has a 
unique inverse (under the condition that <p^'^^ is cyclic) and thus (j)^"^^ is expressed in terms 
of (^("^^ with m < n. 

A simpler iterative scheme can be developed by starting with the (trivial) equation 

4> = 1 + uidi(f) + U2d24> (16) 
and using ([T2|) to eliminate the terms with 82 derivatives, producing 

(p=l + uidicp + ^U2 {di(p - gdfcp - (puicp) . (17) 

Given (j)^"^^ as in (|T^), this equation determines the terms in ^("^+^) with an initial U2 (by 
simply dropping the terms 1 +uidi(j)). There is then a unique cyclically symmetric (1)^'^+^^ 
which contains such terms. 

It may be interesting to note that a purely algebraic equation can be derived for (J). 
We start with (|1^) and its di derivatives 

(p = 1 + UidiCp + U2d2(p 

di(t> = di(/}{0) + uidfcf) + U2d2di(j) (18) 
df(t> = df(t>{0) + uidU + U2d2dl<t). 

Using (|T^) to eliminate the terms with 82 derivatives produces equations linear in 
d^cf) and Of (p. One then takes the system of these three equations and (|T3|) , and eliminates 
these higher derivatives between the equations, producing a single (non-linear) equation 
in terms of (p and di(p. One then repeats the same derivation with the roles of ui and U2 
interchanged, producing an equation in terms of (p and d2(p- Finally, these two equations 
and ([TE| ) can be combined to eliminate dicp and d2(p- We omit the final result due to its 
length. 

These manipulations can be carried through despite the non-commuting nature of wi 
and U2- The final equation however contains the inverses u^^ and which require some 
discussion. First, the combinations wim^^mi and U2u'^^U2 appear. This is incompatible 
with the expansion (|15D and such terms must cancel in the final result. Second, the 
combination U2^<p appears. It is not clear to us that the iterative solution ( |T5D is possible. 

It would be interesting to express the results more explicitly. Of course, we considered 
(^4[) a complete solution of the truncated problem, even though its solution could be made 
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more explicit, because the techniques for solving such equations are so familiar. In practice 
one rarely even writes down its exact solution and instead studies limits of it, or its analytic 
behavior. Perhaps one should regard the solution of (|T3| ) in the same spirit. 

We now turn to study a similar equation for the master field. One such equation 
takes the form 

^ -S[M]=71, + [A,m]. (19) 



To deal with this explicitly we must choose a representation for the master field. Many 
representations exist, in principle related by similarity transformations, and suited for 
different purposes. A very natural one is the 'C-representation,' [^,|^ built from ^^{j), 
the generating functional of connected planar correlation functions, a function of non- 
commuting variables @] The master field Mi is simply 

with the non-commuting derivative as above. If we define a trace as 

tid = di\^^^ = {\d\) (21) 

these operators are a master field. 

The relation between the two generating functionals is 

<^(«) = ^(i), h = uMj)- (22) 



and an equation with the same content as (]T^) was derived in 0] by simply changing 
variables from the Schwinger-Dyson equations. In general it is 

^ S[M]\)=M)- (23) 



Comparing with (|T9|) or the similar equation in this equation is a bit simpler, as one 
bypasses the problem of finding the similarity transformation A (or representing the algebra 
[ni,Mj] = On the other hand, those equations are representation-independent, 

while this one is tied to the C-representation. (The other difference between (|T^ and 
the other equations is that it produces cyclically symmetric equations which sum over all 
variations of Mi in the original path integral, while the others involve a single variation.) 
For our case (^Bp becomes 

diij - cd2^ - gdliP - g{di^f =h (24) 
d2i) - cdiij - gdlij - gid2ipf =j2- (25) 
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We first derive an algebraic equation for V'olii) = '^{31^32 = 0) by a truncation 
analogous to Truncating gives 



aiV'o - c-i/^i - gdiii/jodii/jo) = ji 

(26) 

ditpi - ctp2 - gditpi - gdiipidiipQ - gditjj{0)ditpi = 

while ( p5D gives 

i/^i - cdii^o - gi^2 - gii^i)^ = 0. (27) 

Again we have three equations for three unknowns, and it is easy to see that a quartic 
equation can be obtained for -00- 

In fact this equation can be directly related to (p^. Let K{ji) = l/ji-i/^oOi) and use 
the reciprocity relations ji = R{z), z = K{ji). The algebraic equation F{z,R{z)) = is 
just F{K{ji), ji) = 0. From (p^), it is a quartic equation in terms of K. 

To solve (^3]) more generally, more steps are necessary. Applying 82 to (p^), then 
replacing 91-0 by solving (|2^), we have 

92^1 '0 - gd2dli! - gd2diijdii> - gdii/; [0)828111; 

- - (^2^ - cai^A - j2) + 0(^2^')' = 0. ^^^^ 

g 



One can replace 821P in the above equation by solving the first equation of (p^). To replace 
82811I) and 828^1!^, we apply 81 to the first equation from the right, and use the cyclic 
symmetry of coefficients of "0. Thus, 

821/; =- {8ii; - g8li; - g{8iipf - ji) 
c 

8281^; =^ {8fij - g8li) - g8iij8fij - g8it(;{0)8fij - l) (29) 
828fij =^ {8fij - g8ti; - g8ii^8fi^ - g8fi;{0)8fij - g8ii;{0)8fi;) . 

Substituting these into (|2^), a single equation containing only derivatives with respect to 
ji is obtained. This equation can be solved iteratively in the same way as (^3p. Using (PP]), 



this solution determines the master fields Mi for this model. We intend to study these in 
more detail in subsequent work. 

To summarize, we formulated and solved both large N Schwinger-Dyson equations and 
master field equations in the two matrix model, in terms of functions of non-commuting 
variables. The master field equation had the same basic structure as the Schwinger-Dyson 
equation and the same truncation method worked to solve it. 
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It is interesting to ask if truncation can work in more models. Another model which 
can be solved by the methods presented here has action S{Mi, M2) = tr (M^ + M| — 
CM1M2M1M2) /2, and we will discuss this and the general 'meanders' problem in future 
work. 

There is no reason to think that similar truncations exist for all interesting models 
and it would be very desirable to develop more general techniques to work with this type 
of equation. For physics, it is often better to have qualitative and approximate approaches 
to understanding an equation, and if we can precisely define the equations, we can try to 
study these directions systematically. The results presented here should provide a useful 
test case. 

It is also interesting to see whether one can go beyond the large limit with the 
master field, for example to incorporate considerations in [|13| . 
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